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Abstract 

In this paper, the generalized second law (GSL) of thermodynamics and entropy is revisited in 
the context of cosmological models with bouncing behavior such as chameleon cosmology where 
the boundary of the universe is assumed to be enclosed by the dynamical apparent horizon. From 
a thermodynamic point of view, to link between thermodynamic and geometric parameters in 
cosmological models, we introduce "entropy rate of change multiplied by the temperature" as a 
model independent thermodynamic state parameter together with the well known {r, s} statefinder 
to differentiate the dark energy models. 
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1. INTRODUCTION 



Recent observations of high redshift type la supernovae, the surveys of clusters of galax- 
ies, Sloan digital sky survey ( SDSS) and Chandra X-ray observatory reveal the universe 
accelerating expansion and that the density of matter is very much less than the critical den- 
sity l]. In addition, the observations of Cosmic Microwave Background (CMB) anisotropies 
indicate that the universe is flat and the total energy density is very close to the critical one 

a. 

The observations strongly indicates that the universe presently is dominated by a 
smoothly distributed and slowly varying dark energy (DE) component. A dynamical equa- 
tion of state ( EoS) parameter that is connected directly to the evolution of the energy 
density in the universe and indirectly to the expansion of the Universe can be regarded 
as a suitable parameter to explain the acceleration and the origin of DE m-m. In scalar- 

fin 

tensor theories [5j-|10||, interaction of the scalar field with matter ( for example in chameleon 
cosmology) can be used to interpret the late time acceleration and smoothly varying EoS 
parameter [lll]-fl3|. 

Motivated by the black hole physics, it was realized that there is a profound connection 
between dynamic and thermodynamic of the unvierse (see for example [l^]-j26|). In par- 
ticular, the validity of the GSL j^] which state that entropy of the fluid inside the horizon 
plus the entropy associated with the apparent horizon do not decrease with time, has been 
the subject of many studies. In here, the GSL is discussed for models such as chameleon 
cosmology with bouncing behavior. 

In order to differentiate between cosmological models, a sensitive and robust geometric 
diagnostic for dark energy models is proposed by Sahni et al. 



that makes use of param- 



eter pair {r, s}, the so-called statefinder. It probes the expansion dynamics of the universe 
through higher derivatives of the expansion factor as a natural companion to the deceleration 
parameter. In this work we relate the geometric properties of the cosmological models with 
the thermodynamic one via statefinder parameters. The "entropy rate of change times the 
temperature" is defined in terms of the second derivative of the scale factor of the universe, 
instead of deceleration parameter. This parameter is related to the statefinder parameter, 
and together are adopted to differentiate between dark energy models. 

Section one is devoted to the model independent, thermodynamic formulation of the 
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cosmological models in relation to the dynamical parameters. In section two, we derive 
the field equations for chameleon model and examine its thermodynamic behavior. The 
statefinder parameter from a thermodynamic point of view is presented in section three 
followed by its application in chameleon model and a summary. 



2. GSL AND ENTROPY 

In the following we make two assumptions: 1) in addition to the entropy of the universe 
inside the horizon, an entropy is associated to the apparent /event horizon. 2) with the local 
equilibrium hypothesis, since the temperature of the fluid is almost equal to the apparent 
horizon temperature, the energy would not spontaneously flow between the horizon and the 
fluid inside the horizon. 

According to the GSL in an expanding universe, entropy of the viscous dark energy 
(DE), dark matter (DM) and radiation inside the horizon plus the one associated with 
the apparent horizon do not decrease with time. In general, there are two approaches to 
validate the GSL on apparent /event horizons: 1) by using first law of thermodynamics and 



find entropy relation on the horizons i.e. 



29(], (30] 



T h dS h = -dE h = AnRlHT^^dt = ^R\H{p eff + p eff )dt, (1) 

where k m = (1, —Hr, 0, 0) are the (approximate) Killing vector ( the generators of the hori- 
zon), or the future directed ingoing null vector field and "h" stands for the horizon. 2) in the 
field equations, by employing the horizon entropy and temperature formula on the horizon, 

S h = nRl (2) 

Note that only on the apparent horizon the two approaches are equivalent [? ]. According 
to the recent observational data from type la Supernovae in an accelerating universe, the 



enveloping surface should be the apparent horizon rather than the event one 



311 ] . Therefore, 



re dynamical apparent horizon 



32]. 



from now on, we assume that the universe is enclosed by t 
with the radius given by Rh = ^=f in a flat FRW universe 

By the horizon entropy and temperature given in equations ([2]) and ([3]), the dynamics of 
the entropy on the apparent horizon is 2^], 

S h = 2nR h R h . (4) 



Also, from the Gibbs equation, the entropy of the universe inside the horizon can be related 
to its effective energy density and pressure in the horizon with, 

TdS in = p ef f dV + d(E in ), (5) 

where Si n is the internal entropy within the apparent horizon and p e ff is the effective pressure 
in the model. If there is no energy exchange between outside and inside of the apparent 
horizon, thermal equilibrium realizes that T = Th- Hence, the expression for internal energy 
can be written as E in = p e ffV, with V = ^kR\. From equation dSJ), by using Friedmann 
equation in the cosmological models and doing some algebraic manipulations we find that 
the rate of change of the internal entropy, horizon entropy and total entropy are respectively, 

S in = 12ir 2 R 2 h H(l + Lo eff )(l + 3tu eff ), (6) 
S h = 24ir 2 R 2 h H(l + u eff ), (7) 
S t = 36n 2 R 2 h H(l + uj e f f ) 2 . (8) 

For the rate of change of the internal entropy, equation ([6]), we find that for an expanding 
universe with acceleration, H > 0, and —1 < u) e ff < —1/3 in quintessence era, Si n < 0. On 
the other hand, in phantom era, u e ff < —1, and decelerating universe, u e ff > —1/3, we 
obtain Si n > 0. From equation (J7J), it can be seen that again in an decelerating expanding 
universe and when u e ff > 0, then Sh > 0. Otherwise, Sh < 0. Finally, in equation flSJ), the 
sign of the total rate of change of the entropy, St, for an expanding universe is independent 
of EoS parameter. 

3. CHAMELEON MODEL 

In this section, we consider the chameleon gravity in the presence of matter with the 
action given by, 

S = J - + V{<j>) + /(0)L m ]v/=^b 4 , (9) 

where R is Ricci scalar, G is the newtonian constant gravity and is the chameleon scalar 
field with the potential V(<f>). Unlike the usual Einstein-Hilbert action, the matter La- 
grangian L m is modified as f(<f))L m , where /(</>) is an analytic function of <p. This last term 
in Lagrangian brings about the nonminimal interaction between matter and chameleon field. 
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The variation of the action fl9]) with respect to the metric tensor components in a spatially 
flat FRW cosmology yields the field equations, 

3# 2 =p m + ^0 2 + ^(0), (10) 
2H + 3H 2 = - 1Pm - + V(<j>), (11) 

where we put 8irG = c = h = 1 and p m = p m (a, 0). The dots means derivatives with respect 
to the cosmic time t. The energy density p m is the matter energy density in the universe. 
Also variation of the action Q with respect to the scalar field provides the wave equation 
for the chameleon field as, 

0(0 + 3#0) = -V- (lnX)'p m , (12) 

where we have introduced X(<fr) = /(0) 1-37 . The parameter 7 is the barotropic index of 
the background fluid (DM). From equations (flOl) . (TTTj) and (Tl2|) . one can easily arrive at the 
relation, 

Pm + 3H(l + 1 )p m = (lnX)' Pm <P, (13) 

which for p m (a,<f>) = pom(a)f((j)) readily integrates to yield 

p m (a, 0) = MX^aT 3 ^ = Mr 3 ^a~ 3 ^\ (14) 

with M as a constant of integration. From equations fflOl) and (fTTj) and in comparison with 
the standard friedmann equations we identify p e // and p e ff as, 

Peff = Pm + l^ + V{<P). (15) 
Peff = IPm + ^0 2 - V{4>), (16) 

with the equation of state, p e ff = u e ff Pe ff. 

For the cosmological models including chameleon cosmology with a possibility of bouncing 
behavior, the total entropy rate of change can be negative, thus, GSL is violated when the 
universe is in contraction era. To study the thermodynamic behavior of the chameloen 
model at late time, the structure of the dynamical system is revisited by taking into account 
the following dimensionless variables, 

O — ^ 2 O _ Y_ O c- _ Pmf_ , 17 s 

^ L <j>K.E — QfJ2'> <I>P-E. ~ 3^2' mJ ~ 3f/~2' ^ ' 
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In order to close the system of equations we make the following ansatz: We consider that 
f{4>) = /oexp(5i0) and V(<f) = Voexp(£ 2 0) where 5i and 5 2 are dimensionless constants 
characterizing the slope of potential V(4>) and /(</>). There are no priori physical motivation 
for these choices, so it is only purely phenomenological which leads to the desired behavior 
of phantom crossing. Using equations ( TTUj) -( TTB"j) . the evolution equations of these variables 
become, 

= fi^[V65 2 ^Vl+(l + 7)^/ + 6(Q m /) 2 ], (19) 

<^11 = n m /[-3(i + T ) + 3(i + 7 )n m / + v^«yi(i-37)v^^+6fi 0Jc .j, (20) 

where iV = ln(a). The Friedmann constraint equation fllOp can be written as 

fiu + fi ro / = l, (21) 
where = E + Q<f> PE - In term of the new dynamical variables we also have, 

^ = ^[3(l + 7)n m / + 6n^„]. (22) 

which can be used later in deriving other dynamical parameters. Using the Friedmann 
constraint ( )2~Tj) . equations ( !T8|) - (!20|) reduces to 



dN 



2n feE [3(i - 7 ) + ^(1 - 3 7 )5ifi m / - 6(1 - n ro / - 



-3(l- 7 )>/*V3. (23) 
= fi m /[-^^ v / ^-3(l-fi m /-^,J + ^ 2 ] 

- ^ (1 a (24) 

It is more convenient to work with the equations ( 1231) and ( |24|) than equations ( fT8l) -(l20l). 



With respect to these equations, the total entropy rate of change becomes, 

^ = 87r ( j^) 2 = 2^(3(1 + 7)OrJ + 6^J 2 . (25) 

With numerical computation, in Fig.l, the dynamics of total entropy rate of change and 
EoS parameter with respect to the redshift z for 7 = 0, 1/3 are plotted. 
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Fig. 1: The trajectories of the total entropy rate of change evolution and 
EoS parameter against redshift z for 7 = 0, 1/3. 
ICs: f20 KJ5 |o = — ^P, ^m/lo — ~|- 

From the graphs, at higher redshifts, the EoS parameter is positive and the universe 
decelerates. The S t gradually increases to its local maximum value when u e ff is maximum at 
z ~ 0.7 in a decelerating expansion universe. It then sharply decreases to its local minimum 
value, (St = 0), correspond to the phantom crossing in the past with a fast transition from 
a decelerating expansion to an accelerating one. By entering the phantom era, (w e // < — 1) 5 
S t starts to increase to its local maximum value when EoS parameter reaches its minimum. 
While EoS parameter continue its journey from phantom era to quintessence era passing once 
again the phantom barrier in future, the S t function reaches its local minimum and then 
increases. Since St function vanishes at its local minimums where St = at the phantom 
crossing, the entropy function has points of inflexion at crossing points. 



4. STATEFINDER DIAGNOSIS 

In general, for the statefinder of the cosmological models, we define u = S t T = 87r(-^-) 2 , 
as a thermodynamic parameter, instead of the geometric deceleration parameter. Thus, the 
statefinder parameter {r, s} as the geo-thermodynamics parameter in terms of u become, 

H 3.n. 1/9 „ 2(r — 1) . 

r = Jp-2^' +1 ' S =3(^3- <26) 
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We therefore have in the three cosmological scenarios, {u, r, s} = {187r, 1, 0} as a fixed point 
for SCDM state, {u,r,s} = {2tt, 1,0} as a fixed point for LCDM state, and {u,r, s} = 
{187T, 1,1} as a fixed point for SS state. In terms of the dimensionless variables in the 
chameleon model we obtain, 

r = 3v/lW3 + ^4(1 - 7 )(1 - 3 7 ) + V65 2 n^ E ) + + 7 ) + 1, (27) 

-6£y g - (1 + 7)^/(37 + 6^(1 - 3 7 ) y/H^J - 2y^ 2 y'^WW, 

6[l + (l + 7 )n ro / + 2fi^.J 

We also have tt in terms of these variables as, 



(28) 



27r[3(l + 7 )fi m / + 6^ 



(29) 



Fig.3 shows the statefinder diagrams {u, z}, {u, r} and {u, s} evolutionary trajectories for the 
chameleon model from a thermodynamic point of view. From the graphs, all the trajectories 
for 7 = 0, 1/3 commence evolving from a critical point in the past and end their evolution to 
a second critical point in the future. The current values of the statefinder on the trajectories 
show that the universe currently is between LCDM state in the past and SS state in the 
future. There is thermodynamic stability during the universe expansion so that the universe 
cools down through the expansion. The universe transits from a state in the past to the 
final state in the future while not passing through the above states. 
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Fig. 3: The trajectories of the statefinder {r, u] and {s,u} and u against redshift z, 
for both 7 = 0,1/3. ICs: % K E | = -f , Q m f\ = -§. 

In summary, we realize that the total entropy of the universe in chameleon cosmology 
with bouncing behavior and effective EoS parameter increases/decreases with time in the 
expansion/contraction period thus GS1 is violated in the contraction period. The total 
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entropy rate of change as a model independent parameter is discussed for the model. For 
the chameloen model, the universe starts from a state in the past with St > 0, excluding 
the early early time when a possible bouncing occurred, and following its journey to another 
state in future again with positive S t . Phantom crossings in the past and future correspond 
to the S t \min — when the entropy function has points of inflexion. We then introduce S t T 
to link between geometric and thermodynamic behavior of the cosmological models, instead 
of deceleration parameter. This parameter together with the corresponding statefmder {r, s} 
can be used to differentiate the dark energy models from a thermodynamic point of view. 
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